ACT-09/96 
TAC-96-015 
UCSBTH-96-19 
UM-AC-96-05 
October 1996 
hep-ph/9610405 



BARYOGENESIS DURING 
REHEATING IN NATURAL INFLATION 

AND COMMENTS ON SPONTANEOUS BARYOGENESIS 

Alexandre Dolgov* 
Teoretisk Astrofysik Center 

Juliane Maries Vej 30, DK-2100, Copenhagen, Denmark 

Katherine Freese^ 
Physics Department, University of Michigan, Ann Arbor, MI 48109 

Raghavan Rangarajan^ 
Astroparticle Physics Group, Houston Advanced Research Center 
The Mitchell Campus, The Woodlands, TX 77381 

Mark Srednicki 

Department of Physics, University of California, Santa Barbara, CA 

93106 

ABSTRACT: We calculate the baryon asymmetry created by the decay of a 
pseudo Nambu-Goldstone boson whose interactions violate baryon number con- 
servation. Our results are in disagreement with previous results in the original 
spontaneous baryogenesis models for the asymmetry produced by the decay of an 
oscillating scalar field with B number violating derivative couplings; we find that 
the net baryon number density is proportional to Of, where 9^ is the amplitude 
of the PNGB-field in natural inflation at the onset of reheating. We also discuss 
our disagreement with the interpretation of 9 as an effective chemical potential 
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Section 1: Introduction 

In this paper, 1) we calculate the baryon asymmetry obtained during reheating follow- 
ing natural inflation using an approach that is generally valid for baryogenesis models using 
decaying classical fields. Our results are in disagreement with the the results presented in 
the original spontaneous baryogenesis papers. 2) We discuss an objection to the effective 
chemical potential interpretation used in models of spontaneous baryogenesis. 

In natural inflation the role of the inflaton is played by a pseudo Nambu-Goldstone 
boson, hereafter referred to as 6>, with a potential of the form [1] 

V(9) = A 4 (l-cos6>). (1.1) 

This model was proposed to "naturally" provide the flat potential required for inflation 
to work [2], [3]. Here 9 = $//, where $ is a complex scalar field and / is the scale at 
which a global symmetry is spontaneously broken; soft explicit symmetry breaking takes 
place at a lower scale A. From eq. (1.1) one can see that the height of the potential is 2 A 4 
while the width is /. Since the scales of spontaneous and explicit symmetry breaking can 
"naturally" be separated by several orders of magnitude, one can obtain A < 10 -3 / as 
required for successful inflation [4]. 

In ref. [14] an extensive study of the conditions under which the 9 field can drive 
inflation has been obtained. After the period of inflation, the energy density of the 9 field 
is converted to radiation during reheating through its decay to other forms of matter as 
it oscillates in its potential. Below we shall assume that 9 is coupled only to fermions. 
We treat 9 as a classical scalar field coupled to quantized fermion fields Q and L via an 
interaction term of the form QLe^ + LQe~ , where Q carries baryon number but L does 
not. We show that the decay of 9 gives rise to a net baryon number density (n& — n^) 
proportional to 6> 3 , where 6>j is the value of the 9 field at the onset of reheating. 

Our result disagrees with the calculation in the original spontaneous baryogenesis pa- 
pers [15] where it was argued that the asymmetry is proportional to 6>j to the first power, 
independent of the details of the baryon number violating couplings of the 9 field. Specif- 
ically, in previous work, Cohen and Kaplan [15] considered any theory in which a scalar 
field is derivatively coupled to the baryon current with a term in the interaction La- 
grangian of the form C- m i oc 8^9 J^ 1 , and derived an expression for the baryon asymmetry 
produced by the decay of the scalar field as it oscillates about its minimum. The pseudo 
Nambu-Goldstone boson in natural inflation can serve as an example of such a scalar field. 
Cohen and Kaplan obtained |ng| = r/ 2 |#|, where Y is the decay rate of the 9 field and tlq 
is the net baryon number density. This gives 

|Ang| = Tf 2 \A9\ . (1.2) 
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Below we discuss our concerns with this conclusion and present calculations for the specific 
case of eq. (1.1); our results disagree with eq. (1.2). We also comment on our objections 
to interpreting 9 as a chemical potential when 9 is small, as was done in ref. [15]; we 
argue that a Lagrangian term 9 J° does not appear in the Hamiltonian, and therefore it is 
incorrect to identify 9 with an effective chemical potential for baryon number. 

The framework of this paper is as follows. In Section 2, we write down the Lagrangian 
density for the inflaton field and present the equation of motion for 9 as it oscillates during 
the reheating phase, as derived in ref. [16]. In Section 3 we discuss our concerns with 
eq. (1.2) as obtained in ref. [15] (these concerns were raised in an earlier paper [16] by 
two of the authors [Dolgov and Freese]). We then proceed to calculate the total baryon 
number and antibaryon number produced during the decay of the 9 field, and find a 
net baryon number density (n^ — n^) proportional to 9f. We also show that the energy 
density of the produced particles is equal to the initial energy density of the 9 field as 
a a check on our calculation. In Section 4, we discuss how constraints on parameters in 
natural inflation obtained in ref. [14] affect the quantitative results for baryogenesis. We 
also discuss our objections to the thermodynamic generation of the baryon asymmetry 
via an effective chemical potential interpretation in models of spontaneous baryogenesis. 
Finally we summarize our results. In the Appendices we provide details of the calculations 
outlined in the main body of the paper. In particular, in Appendices A and B, we include 
derivations of the number density of particles produced by the decay of a classical scalar 
field; the number density of particles produced is proportional to the integral over momenta 
of the one pair production amplitude. 

Section 2: The Model 

As in ref. [16] we consider a simple model involving a complex scalar field $ and 
fermion fields Q and L with the Lagrangian density 1 

£ = -d^d^-V(^^)+iQ^d^Q+iL^d^L-m Q QQ-m L LL + (g^QL+h.c.) . (2.1) 

Note that, despite their names, Q and L cannot be actual quarks and leptons, since the 
interaction term does not conserve color. They could, however, represent heavy fermions 
with other interactions with the fields of the Standard Model which fix the assignments of 
global charges. In particular, we shall assume that the field Q carries baryon number while 
the field L does not. The U(l) symmetry that corresponds to baryon number is therefore 
identified as 

$^e m $, Q^e ia Q, L^L. (2.2) 

We use a metric (-1,1,1,1). 
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We assume that this global symmetry is spontaneously broken at an energy scale / via a 
potential of the form 

V(\$\) = A($*$-/ 2 /2) 2 - (2.3) 
The resulting scalar field vacuum expectation value (VEV) is (<&) = fe^lf 



Below the scale /, we can neglect the radial mode of $ since it is so massive that it 
is frozen out; w ra( jj a j = A 1 / 2 /. The remaining light degree of freedom is 0, the Goldstone 
boson of the spontaneously broken U(l). For simplicity of notation we introduce the 
dimensionless angular field 9 = <fi/ f. We then obtain an effective Lagrangian density for 
6, Q, and L of the form 

f 2 

C eS = -—dpfid^d + xQ^d^Q + iL^dfiL - mnQQ - m^LL + {-^=fQLe %6 + h.c.) . (2.4) 
2 v2 

The global symmetry is now realized in the Goldstone mode: £ c g is invariant under 

Q^e ia Q, L^L, 9^9 + a. (2.5) 

With a rotation of the form in eq. (2.5) with a = — 6>, the Lagrangian can alternatively be 
written 

f 2 - - Q 

C eff = - J —d^9 + iQ^d l ,Q + iL^d^L-mQQQ-m L LL+(-^fQL + h.c.) + d fl 9J' 1 , 

(2.6) 

where the fermion current derives from the U(l) symmetry; here, = Q'j^Q. 

We now assume that the symmetry (2.2) is also subject to a small explicit breaking 
which gives rise to a potential as in eq. (1.1) and which provides a nonzero mass for the field 
9. This explicit symmetry breaking could come from Planck scale physics. Alternatively, 
one can imagine a scenario similar to that involving the QCD axion where, at energy scales 
of the order of Aqc D , instanton effects create the fermion condensate (ipij>) ~ AqcD' giving 
rise to a mass term for the axion. Note that for the natural inflation model, the required 
mass scales are much higher than for the QCD axion. The width of the potential must be 
roughly the Planck mass in order to achieve enough e-foldings of inflation, and the height 
of the potential must be roughly A 4 ~ [10 16 GeV] 4 in order for density perturbations 
appropriate for structure formation to be produced (see the Discussion section at the end 
of the paper for more detail). Consequently the scale at which the relevant gauge group 
(not QCD) must become strong is roughly the GUT scale. These and other mechanisms 
such as those found in technicolor and schizon models for generating a potential for pseudo 
Nambu- Goldstone bosons are discussed in ref. [4,14]. 



Initially, as the 9 field rolls down towards the minimum of its potential, its potential 
energy drives inflation. Let 6^ be the value of the 9 field at the beginning of the reheating 
epoch, after inflationary expansion has ended. (We shall ignore spatial variations in the 9 
field.) During the reheating epoch the 9 field oscillates about the minimum of its potential. 
While 9 oscillates it decays to the fields Q and L. The interactions of the fermionic fields 
create a thermal bath thereby reheating the universe. Note that we must take g <C 1 so 
that fermion masses generated for the fermions from the Yukawa coupling, ~ gf, are 
small enough that the fermions can in fact be produced by decays of the pseudo Nambu- 
Goldstone bosons. See ref. [10] in Dolgov and Freese [16] for further discussion of this 
point. 

The equation of motion for the 9 field with the back reaction of the produced fermions 
was rigorously derived in the one loop approximation in ref. [16]. For small deviations 
of 9 from the equilibrium the potential can be approximated as V(9) = ^m^/ 2 6* 2 and 
the equation of motion during the oscillating phase can be effectively written in the well 
known form 

9 + m\9 + T9 = , (2.7) 

where is the renormalised 9 mass defined as lim^^oo m? R [l + £^ ^°s(^ w / m R)] = m2 ; 
where m is the bare mass of the 9 field, and V = g 2 mn/8iT. (Our expressions above differ 
by a factor of 2 from those in ref. [16] because a factor of l/y/2 was dropped from eq. (2.5) 
in ref. [16].) The solution to this equation is 

9(t) = 9 i e~ Tt /' 2 cos(m R t) . (2.8) 

where we have assumed that the initial velocity of the 9 field is negligible and have therefore 
set an arbitrary phase in the cosine to zero. The results obtained below can be easily 
generalized for arbitrary initial conditions. The above solution was derived assuming uiq = 
= 0. However, it can be shown that non-zero values of uiq and will not change 
the solution for 9 significantly as long as mg, <C mji, which we shall assume below. 

Section 3: Baryogenesis 

Previous Calculations and Concerns: In previous work, Cohen and Kaplan [15] consid- 
ered any theory in which a scalar field is derivatively coupled to the baryon current with 
a term in the interaction Lagrangian of the form oc d^9J^, and derived an expres- 
sion for the baryon asymmetry produced by the decay of the scalar field as it oscillates 
about its minimum. From eq. (2.6) one can see that our pseudo Nambu-Goldstone boson 
is an example of such a scalar field as it has the appropriate coupling. Cohen and Kaplan 
obtained |ng| = r/ 2 |#|, where rig is the net baryon number density. This gives 

|Ang| = r/ 2 |A6>| . (3.1) 
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In a previous paper [16] by two of the authors [Dolgov and Freese], several concerns with 
this interpretation were raised. We will outline two of these concerns again here, and then 
proceed with a direct calculation of the baryon asymmetry. Our results will disagree with 
eq. (3.1). 

One concern is as follows: in making the identification |ng| = r/ 2 |#|, one is comparing 
an operator equation, namely, the Euler-Lagrange equation 9 + w?9 = ng/ / 2 , with an 
equation of the form of eq. (2.7) which is obtained after vacuum averaging. In ref. [16] the 
average value (ng) was found to be not just —Tf 2 9 but a more complicated expression 
(eq. (3.3) in ref. [16]). 

A second concern is with regard to energy conservation. The initial energy density of 
the field 9 which creates the baryons and antibaryons is po(ti) ~ / 2 m 2 6> 2 . At the end some 
of this energy density has been converted to baryons and antibaryons, with energy density 
p(tf) > tiqEq where Eq ~ m is the characteristic energy of the produced fermions 
(note that rig refers to the difference between baryon and antibaryon number densities 
and not to the total number density of produced particles). Clearly it must be true that 
p(tf) < Pe(ti). If we were to use eq. (3.1) we would see that this requires V < A9m. From 
the definition of V we see that this is satisfied for small values of coupling constant g as 
long as A9 is not too small; for small values of A6>, this relationship can never be satisfied. 

New Calculations and Results: We now proceed to calculate the net baryon number 
density of the particles produced during reheating. We perform an explicit calculation 
and find a different result from eq. (3.1). The 9 field decays to either QL pairs or QL 
pairs. (The Q and L fields are not the mass eigenstates. Later in this section we consider 
effects of oscillations between Q and L fields.) As mentioned earlier, we treat the 9 field 
classically, Q and L are quantum fields and Q carries baryon number. For now we ignore 
any dilution of the baryon number density due to the expansion of the universe. 

As shown in Appendix A with the Bogolyubov transformation method [17], the aver- 
age number density n of particle antiparticle pairs produced by decay of a homogeneous 
classical scalar field, to lowest order in perturbation theory, is given by 

n = ^J2 fdPldP2\A\ 2 , (3.2) 

where A is the one pair production amplitude, subscripts 1 and 2 refer to the final particles 
produced and dp = d^p/[(2n)^2p°]. Eq. (3.2) can also be obtained using the method 
presented in Sec. 4-1-1 of ref. [18], as discussed in Appendix B. 

Thus, to lowest order in perturbation theory, the average number density of QL pairs 
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produced during reheating in our model is given by 2 

n(Q,L) = l [dpdq\{Q(p,s Q ),l(q,s L )\0)f 

We take 



(3.3) 



(3.4) 



and a similar expression for L. Here {bLbfi} = {di,d s '}} = (2n) 3 2k° 5 3 (k - k')<5, 



ss' 



Standard algebra gives 
n(Q, L)= y Yl J d P d( l 

SQ,Sl 

9 2 f 2 



■ 9 



(Q(p,s Q ),L(q,s L )\i^=Jd*xQ(x)L(x)e ie W\0) 



2V 



dp dq 



/oo 
dt j2wt+i6{t) 
-co 



Tr[(-ff + mQ){-4 - m L )} 

(3.5) 

where 2uj = p + y . We obtain a similar expression for n(L,Q) with 9{t) replaced by 
—9(t). We set the baryon number density to be equal to n(Q,L) and the antibaryon 
number density to be equal to n(L, Q). Then we have 



2/2 



n 



6,6 



2tt s 



eft e 2^ e ±<«(*) 



■oo 



(3.6) 



where the + sign in the exponent refers to baryon number and the — sign to antibaryon 
number. To carry out the integral over time we expand e 1 ^ as 



l + i6-6 2 /2, 



valid for small 6, and use 



0(f) 



for t < 



9ie~ Tt / 2 cos(m R t) for t > 



(3.7) 



(3.8) 



We also use a convergence factor at early times to regularize the integral. We will examine 
a series of possible terms to find the first nonzero contribution in perturbation theory. 
The lowest order term comes from using e l ® = 1 from eq. (3.7) in eq. (3.6) and gives 
/ dte 2lLut oc 5{2u) = since we can not have to = for particle production. The next term 
in the expansion, the 9 term in (3.7), when squared gives the same contribution to and 



2 Throughout the paper, a state (A(p\, s\), B(p2, S2)\ corresponds to a final state with an 
A particle of momentum p\ and spin s\ and an anti-B particle with momentum p2 and 
spin S2- 



to rift. In order to obtain an asymmetry one must consider cross terms. The lowest order 
cross term that gives a nonzero contribution to the baryon asymmetry is 



2x2 



nz = 2 x 



2^ 



0(2u) 6 2 (2oo) 



2i 



+ h.c. 



(3.9) 



where h.c. refers to hermitian conjugate, 

9{2u) = J 

and 



dte 2lujt 6(t) 



DO 



9 2 (2u) = I dte 2lujt 9 2 (t). 



(3.10a) 



(3.106) 



— oo 



The factor of 2 in eq. (3.9) arises from the fact that the cross terms in n b and n b terms 
are the same up to a minus sign. One can see from the form of eq. (3.9) that we expect 
the asymmetry to be proportional to 6> 3 . The details of this calculation are outlined in 
Appendix C, and the results are presented here. 



We obtain 



1 



2 Q 2 

i 



9 ... j-2n3 



32tt 



Therefore, 



n B = n b 



n- b = 



9 

16tv 



2n3 



(3.11) 
(3.12) 

(3.13) 



= 7?rn 



We notice that the net baryon number density is proportional to 9f. This disagrees 
with the calculation in ref. [15] which gives an asymmetry proportional to 0j. We also note 
that the number density of pairs of particles n^ + n^ is equal to ■^mjif 2 9 2 . Since the energy 
per pair of particles is thr, the energy density in the produced particles is ^m^/ Of, which 
agrees with the initial energy density of the 6 field. We have also done the calculation of 

Pfinal = ^ E fdpdq{p° + q°)[\ (Q(p, sq), L(q, s L ) |0) | 2 + | (L(q, s L ) , Q(p, s Q )\0)\ 2 } 

(3-14) 

and have verified that we obtain \^ 2 jif 2 d 2 - 

Mass Mixing: In many cases eq. (3.13) is not yet the complete story because of mass 
mixing. As we mentioned earlier the Q and L fields are not mass eigenstates. Therefore a 
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particle which is produced as a Q may later rotate into an L. This effect must be taken 
into account. Eq. (3.13) is completely correct for the case where the fermions Q and L 
are converted immediately to regular quarks q and leptons / as soon as they are produced 
(assuming that the temperature is low enough that the q and / cannot convert back into 
Q and L). In that case, there is no opportunity for mixing to take place, e.g., there is 
no opportunity for Q to convert to an L. On the other hand, if Q and L do not decay 
immediately into stable lighter mass particles with appropriate quark quantum numbers, 
they may have the chance to mix into one another. One can calculate the effects of mixing 
in either the Q, L basis or in the basis of mass eigenstates; below we will do both. 

The mass matrix in the (Q, L) basis is 



The mass eigenstates are 



-gf/y/2 m L J 



V>1 = 7= a and ip 2 = (3.16) 

VI + e 2 VI + e 2 



with masses mq — gf/(y/2e) and m^ + gf/(y/2e) respectively, where e = \[2gf /(Am + 
a/ (Am) 2 + 2g' 2 f' 2 ) and Am = mq — m^. Note that Am = corresponds to e = 1. 

In the 1^1,1^2 basis, one can now calculate the baryon asymmetry as a sum of terms, 
each of which is a product of a number density of produced particle/antiparticle pairs 
times the (time averaged) quark content of the pair, 

n B = n(V>i,^ 2 )KWl>| 2 + ^2,V'l)KQI^2)| 2 -n(^i,^2)KQI^2)| 2 -^2,V'l)KQI^l)| 2 - 

(3.17) 

Here n(ipi,ip2) an d ^(V'2)V , l) are the number densities of ifti and ^2 pairs and ^2 an d "4>\ 
pairs respectively; and |(<5|^i)| 2 is the probability that a particle which is produced as a 
ipi (where i = 1,2) is measured as a Q. Hence, for example, the first term is the product 
of the number density of ^2 pairs produced times the quark content of tpi. 

Note that we are here computing a time averaged baryon asymmetry; actually the value 
of the baryon asymmetry oscillates in time, as discussed in Appendix D. From eq. (3.16) 
we see that the probability that V>l,2 is measured as a Q is 

KWPHMI 2 ^ ( 3 - 18 ) 



and 



\(Qm\ 2 = \(Q\M= ^ ■ (3-19) 



As in eq. (3.2), the number densities of particle/antiparticle pairs are obtained by 
squaring the production amplitudes for the pairs, 

n fj = 7 E [J ^ dkj \Af > (3-20) 

where i and j are either 1 or 2. The amplitude for production of a ipiipj pair is 

A fj = J d A x(^=fe id QL + h.c.)\0). (3.21) 
Using eqs. (3.18), (3.19), and (3.20), we can write eq. (3.17) as 

± Wdfc dfo r—yj UA^l 2 - \A 2l \ 2 ] . (3.22) 

Using 

® L = ^ YT? ) $ 2 ^ 1 ~ e ^ 2 ^ 2 + ^i^ 1 _ ^V'l^] (3.23) 
and its hermitian conjugate, we calculate the relevant production amplitudes: 

A 1 2 = (^iMi J d A x(^=fe id QL + ^=fe- id LQ)\0} (3.24a) 



n B = y 



to find 



^2 = 1-2=/^—^ (Vi,^ 2 | | A(M^-^M)|0>. (3.246) 

Now the two matrix elements in eq. (3.24b) are similar to the ones we calculated in eq. (3.5), 
with QL replaced by ^Y^i- Hence, we have 



A 12 = [ ^ ) ( A LQ - ^A QL ) ■ (3-25) 



Similarly, 

= h.rJA^] = I 

1 + e 



A 2l = h.c.[Ad = [ ) (~e 2 A LQ + A QL ) . (3.26) 



Thus eq. (3.22) becomes 

2 



nB = 1 IT? ) ( ^ I dkQ dk ~ L |A ° Z ' 2 ~ ^ I dkL dkQ 

/ \sq,s l J S L ,SQ J 



9 \ 2 
1 -e 2 



(3.27) 



- — 2 j x our previous answer 
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Thus we find that 

If mq = rriL, e = 1 and the asymmetry vanishes because in this case the net baryon number 
of a (-01 , V'2) P arr or a (V^jV'i) P a i r is and thus no baryon asymmetry is produced. 

Another derivation of eq. (3.28) is given in Appendix D. In the preceding paragraphs 
we considered particle production and mixing in the mass eigenstate V^l ? ^2 basis. In 
Appendix D we work in the Q, L basis. We find the oscillations of the baryon asymmetry 
with time, and obtain the same expression as in eq. (3.28) for the time averaged baryon 
asymmetry. 

Thermalization: After the 9 field has decayed into ifri and ^2 particles, thermal equi- 
librium can be established if these particles have other interactions with each other and 
with other particles. As long as one introduces interactions such as xV'lV'l an d X'4 ) 2' l l ) 2 as 
a part of a realistic model, the number of ipi — tf)i particles and of 1^2 ~ V>2 particles does 
not change, thereby preserving the baryon asymmetry. (Interactions such as xV'iV'2 + h.c. 
would, however, destroy the baryon asymmetry.) The fields ipi and ip2 will annihilate 
or decay to lighter particles which will thermalize. If these interactions preserve the net 
baryon number, then the asymmetry will survive. 

Quantitative Results: So far we have not included the effects of the expansion of the 
universe. For baryon number created when H < T, we may neglect the expansion and 
directly use the results obtained above in eq. (3.28) with 9i replaced with the value of 9 at 
H = r. Since the 9 field dominates the cosmic energy density, the condition H = V fixes 
the amplitude of 9 at that moment to be 

01 = y^/4^(Tm Pl /fm R ) « 0.02g 2 m Pl /f « 1. (3.29) 

In the early stages of reheating with 9 > Q\, expansion of the universe must be taken into 
account. 

The decay of the 0-field produces relativistic ipi^ and ipi^ with energies uj ps rnji/2. 
This state is far from thermal equilibrium (the temperature of the thermalized plasma 
in eq. (3.30) below may be smaller than the ip masses). The rate of thermalization de- 
pends upon the interaction strength of the fermions created in the 9 decay. It is typically 
higher than the decay rate because g 1 to ensure reasonable fermion masses. Ther- 
malization could occur either through annihilation of ip\ and ty\ or ^2 and ^2 i n t° light 
particles or through their decays and subsequent elastic scattering. Assuming that these 
processes are fast we can roughly estimate the reheat temperature in the instantaneous 
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decay approximation, p ra( j = p@(t = T 1 ), as 

T reh = (90/8 7 r 3 ^) 1/4 v / Tmp/ » 0.15gg~ 1/A Ay/m P i/f (3.30) 

where we have taken m^> = A 2 //. 

The entropy density after thermalization is given by s = An^g*T^ ^/90. It is conserved 
in the comoving volume if the expansion of the universe is adiabatic, in particular in the 
absence of first order phase transitions as the universe cools. Baryonic charge density is 
also assumed to be conserved inside a comoving volume during and after thermalization 
and so the baryon-to-entropy ratio rig / s remains constant in the course of expansion. 

First we find the baryon asymmetry produced after H < V so that expansion may be 
neglected (subscript 1 refers to this case). Using eqs. (3.28), (3.29) and (3.30) we find 

In the models studied in ref. [14] (f/mpi) = 0(1) and f/A = 10 6 - 10 3 , so to get a 
reasonable baryon asymmetry we need a rather large coupling, g > 10~ 2 (for e < 1). 

In fact the asymmetry should be noticeably larger than that given by eq. (3.31). The 
result that we got above refers to the case when H < V but the process of particle pro- 
duction starts much earlier when H m wir and the inflaton field begins to oscillate around 
the bottom of the potential. The net baryon number density produced while H > V is 
again proportional to # 3 , as it is associated with the interference between the 9 and the 
# 2 terms in | J dte^ iujt (l + i6 — # 2 /2)| 2 in eq. (3.6). The generation of the asymmetry 
is more efficient at early times (H > T) since the amplitude of the 6>-field, which goes 
down with the scale factor as i? -3 / 2 , is larger. However when H > T one must include 
the effects of the expansion of the universe on the production of the baryon asymmetry. 
This makes the exact calculations considerably more complicated. Still we can roughly 
estimate the asymmetry in the following way. The difference between the production of 
particles and antiparticles is most profound at early times, At a ~ I/wir, when 9 is large. 
The total number of particles produced in time At a is proportional to TAt a nQ and, as 
we mention above, the baryon number asymmetry must vary as 9 . Therefore, a reason- 
able estimate of the net baryon number density created while H > V is rig ~ Ff 2 9f. 
Between the time of peak production of baryon asymmetry at t a ~ and the peak 

entropy production at ~ 1/T we will take the baryon asymmetry to be diluted by a 
factor of (Ra/Rb) 3 ~ {ta/hf ~ (r/m^) 2 due to the expansion of the universe, where 
we have taken the universe to behave as matter dominated with R oc in the usual 
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fashion during reheating. Thus the baryon-to-entropy ratio at time % and afterwards is 
(^) 2 ~ T/ 2 6' J 3 (r/m^) 2 /s. The calculation of the entropy density is exactly the same as 
described above eq. (3.31), while the baryonic charge density is larger than the H <T case 
by a factor of (6i/6i) 3 (T/m R ) 2 = 9i/9\ = m R /Y = 8n/g 2 > 1. Consequently, we get that 
the total baryon asymmetry of the universe is approximately equal to 

e^le>.— 3 |!(?) 3/2 f(^) 2 <"*> 

Here subscript 2 refers to the case where expansion has been included. Henceforth we use 
eq. (3.32) as our estimate of the baryon asymmetry produced. 

Section 4: Discussion 

In ref. [14], the authors obtain constraints on the parameters A and /. The stipulation 
that a large fraction of the universe after inflation have inflated by at least 60 e-foldings 
gives / > 0.06Mp;. A stronger constraint can be obtained by requiring the formation of 
galaxies to take place early enough in the history of the universe; in this way one obtains 
/ > 0.3Mpi. A constraint on A is derived by using COBE data on the density fluctuation 
amplitude and is plotted in fig. 1 of ref. [14]; the upper bound on A thus obtained ranges 
from 10 13 GeV to 10 16 GeV for / between 0.3 M P1 and 1.2% If one desires the density 
fluctuations from inflation to be responsible for the large scale structure of our universe 
and hence for the COBE anisotropy, then A must be equal to the above values rather than 
simply being bounded by these numbers. 

If the baryon asymmetry produced above is accompanied by an equal lepton asymme- 
try, so that B — L = 0, it will be wiped out by baryon number violating sphaleron processes 
unless the reheat temperature is below 100 GeV. The low reheat temperature condition 
may be a desirable feature of our model as many inflation models have difficulty creating 
a high reheat temperature. Furthermore, we shall require that T re h > 10 MeV so that we 
reproduce standard nucleosynthesis. If, in addition, one requires the density fluctuations 
from inflation to serve as the explanation for the COBE data rather than merely being 
bounded by it, then A is determined as a function of / as described in the previous para- 
graph; then the combination of these constraints implies that 10~ 14 < g < 10~ 10 for A 
and / equal to 10 13 GeV and 0.3Mp^ respectively, and the asymmetry generated by the 
mechanism considered above is by far below the necessary observed value. However if A 
is merely bounded by COBE measurements (density fluctuations must then be generated 
some other way than by the inflation), then g can be much larger as can the baryon asym- 
metry. Alternatively if a nonzero (B — L) is generated, for example, if the L fields carry 
no lepton number, then it is not destroyed by the electroweak processes and the coupling 
constant g need not be so small. 
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In our perturbative calculations of the number of pairs of particles produced we have 
assumed that the masses of the fermions are smaller than the mass m^> of the theta-field 
and that gf < rriQ ^; otherwise the perturbative approach is not applicable. This implies 
that gf < mji = A 2 // or g < (A//) 2 . In this case, the baryon asymmetry is rather small 
as (^f) 2 < 10 _3 (A//) 5 (mp///) 1 ' 5 < 1CT 18 (in obtaining this limit we have included the 
simultaneous constraint on A and / from density fluctuation constraints in ref. [14]). If, 
however, 9 is not the inflaton field, as in the original version of the spontaneous baryogenesis 
scenario [15], then the parameters A and / do not necessarily satisfy the above bounds 
and the asymmetry may be quite large, especially if / <C trtpi- In such a case, one would 
have to redo the calculation of the entropy if 9 does not dominate the energy density of 
the universe when it decays. A period of inflation prior to the decay of the PNGB would 
also be required so that 9 and, consequently, the baryon asymmetry have the same sign 
within present-day domains of sizes 10 Mpc or greater. (Existing data do not rule out a 
matter symmetric universe with domains of matter and antimatter on scales of 10 Mpc or 
more [19].) 

An interesting possibility is that the mass of fermions is not below mjj and the pertur- 
bative approach is not applicable. The non-perturbative calculations in this case are more 
complicated and will be presented elsewhere. 

We would also like to point out an objection to the mechanism of creating the baryon 
asymmetry ther mo dynamically, via an effective chemical potential interpretation, as first 
proposed in ref. [15] and later applied to spontaneous baryogenesis models at the elec- 
troweak phase transition [20]. The approach in ref. [15] is to identify 9 in the term 
d^JI 1 = -9J° = -9n B in eq. (2.6) with an effective chemical potential that biases the B 
violating interactions in the universe to create more baryons than antibaryons. If this is 
the case then, as is argued in ref. [15], the net baryon number density in the thermal bath 
is given by the expression: rig ~ #T 2 . However this can not be true because it contra- 
dicts the equation of motion of the Goldstone field: <9 2 6> = —d^J^/f 2 . Assuming spatial 
homogeneity, this equation gives rig ~ #/ 2 . (A similar contradiction is obtained using 
the equation of motion for a PNGB-field in an expanding universe.) This contradiction 
is resolved because 9 does not shift energies of baryons and antibaryons and cannot be 
identified with a chemical potential. While the term d^J^ exists in the Lagrangian density 
in eq. (2.6), it does not in the Hamiltonian density 

H = ^fii - L eff (Mi) . (4- 1 ) 
dfa 

where fa represents all the fields in the Lagrangian density [21]. This is similar to the 
interaction of a charged particle with a magnetic field, where the energy of the particle is 
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not changed due to the action of the field as the force is proportional to the velocity and 
orthogonal to it. Thus the term —9n B does not appear in the Hamiltonian density and 9 
can not be interpreted as an effective chemical potential. 

As mentioned above, the approach of ref. [15] has been applied to create the baryon 
asymmetry in spontaneous baryogenesis models at the electroweak phase transition [20]. 
The role of the 9 field is associated with the Higgs field in electroweak baryogenesis. Since 
we feel that the identification of 9 as an effective chemical potential is incorrect, these 
models too are subject to the same criticism. 

In conclusion, we have calculated the baryon asymmetry created by a pseudo Nambu- 
Goldstone boson with baryon number violating couplings in the context of natural inflation. 
We have obtained a general result for the baryon asymmetry created by the decay of an 
oscillating scalar field with baryon number violating couplings and demonstrated explicitly 
that the asymmetry is not proportional to 9{ to the first power as claimed in earlier 
work. We have also discussed our objection to the thermodynamical generation of the 
baryon asymmetry using an effective chemical potential approach in models of spontaneous 
baryogenesis. 
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Appendix A: Number Density of Produced Particles in Terms of One Pair 
Production Amplitude 

Here we use the Bogolyubov transformation method to obtain eq. (3.2). We show 
that in the lowest order of perturbation theory, the average number density of particle 
antiparticle pairs produced by decay of the initial scalar field is given by 

n = Yl I dpi dp2 I" 4 ' 2 ' 

Sl,S 2 

where A is the one pair production amplitude and subscripts 1 and 2 refer to the final 
particle and antiparticle produced. For simplicity we will work with scalar fields here; the 
generalization to production of fermions is similar and has been performed in ref. [22]. 

We begin with a classical scalar field <f>(t) coupled to a quantum complex scalar x- 

L int = 9<t>(t)x*X- (A.1) 
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At early times t — > — oo, we take 0(£) = so that x is expanded in terms of creation and 
annihilation operators, 



X 



/aKexp(-^ + lk .x) +4 texpM- 8k .x)], (A2) 



where uj = Vk 2 + m 2 . Here the commutators are [a kl ,a k2 ] = (27r)^2fcj5^)(ki — k2) and 
a similar relation holds for the antiparticle creation and annihilation operators 6 k . Then, 
at later times, <f>(t) 7^ and eq. (A. 2) is replaced by 

X = Jdk [a k f k (t) exp(zk • x) + exp(-zk • x)] , (A3) 

with equation of motion 

[d? + k 2 + m 2 -g<P(t)]f k (t)=0. (A4) 

The subscript on f k , and on a k and j3 k below, refers to |k| and not to the momentum 
four vector. For continuity at early times f k (t — > —00) = exp(— iut). We also assume that 
4>{i) — > for t — > 00. Then we have 

/ fc (t - +00) -> ajfee"^ + fte* 1 * , (A.S) 

so that x(t) evolves as 

x(t — > +00) = J dk [exp(— iut + ik ■ x)(a k a\ l + /3/|!&L k ) 

+ exp(za;t — ik • x) (a^ + Pk a -k)] ■ 
One can define new creation and annihilation operators, for particles: 



(A6) 



a k = a k a k + Pl b -k ' (A7a) 



and for antiparticles: 



6 k = a k b k + p* k al k . (A.7b) 
Then the operator of final particle number is given by iV k = aj c a k /[2/c°y]. 
The number of particles in the final state of momentum k is given by 

AT k = (0|iV k |0> = |/3 fc | 2 . (A8) 

Thus the total number density of produced particles is 
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This result, obtained by the method of Bogolyubov coefficients, can be found in refs. [17,23] . 

Now we shall calculate (3^ in perturbation theory. Expanding / = /q + /i, we have 
/o = exp(— iujt) and the equation of motion (A. 4) becomes 

(df + k 2 + m 2 )fi = g<f>(t) exp(-iut) . (A10) 



Using the Green's function method we find 



du' cj){u)' — uj) 



W—'J^ J-e-rf ''*"- (A11) 

Taking the residue at the pole uj' = —Vk' 2 + m' 2 = —uj, we find the coefficient of exp(+iujt) 
to be, 

(3 k = tg[4>(2uj)]*/2uj. (A12)) 

Now, for comparison, let us calculate the field theory amplitude with the interaction 
Lagrangian given by eq. (A.l), 

A = (k h k 2 \i j d A x <70(t)x*x|O> • (A13) 

Perturbatively the matrix element is easy to calculate using eq. (A. 2), and we find 

A = z#(27r) 3 5 3 (ki + k 2 ) J dt <p(t) exp[z(u;i + uj 2 )t] , (A14) 

so that 

\A 2 \ = g 2 V(2nf5^\k 1 + k 2 )\4>(uJ 1 + uj 2 )\ 2 • (A15) 
Now if we integrate over dk\ dk 2 , we find that 

1 f d 3 k 2 \4>(2uj)\ 2 , A . 

This is exactly eq. (A. 9) with (3^ given by eq. (A. 12). Thus we have shown that the number 
density of produced particles is given by the integral of the one pair production amplitude 
squared. 

Appendix B: Second Derivation of Number Density of Produced Particles 
in Terms of One Pair Production Amplitude 

Eq. (3.2) can also be obtained using the method presented in Sec. 4-1-1 of ref. [18]. 
(We have ignored the higher order vacuum graphs that give the exponential factor exp(— n) 
in eq. (4-23) of ref. [18].) We have verified that we obtain the Poisson distribution for 
the number of (Q,L) pairs and (Q,L) pairs as in ref. [18]. Indeed the derivation of the 
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Poisson distribution can be done exactly along the same lines as in ref. [18]. The only 
difference is that in the example considered in this book the matrix element describes 
the production of a single photon by an external current while in our case it gives the 
amplitude for production of a pair of particles. For the multiparticle production amplitude 
this gives rise to a different normalization, namely, in the case of the production of n 
photons the amplitude contains the factor 1 / \fn). connected with identical photons while 
for the case of production of n pairs of QL (or charge conjugate) the amplitude contains 
1/nl. In the case of photons the multiparticle amplitude squared contains the following n- 
dependent factors: |A^| 2 ~ |(n!)(l/n!)(l/\/n7)| 2 ~ 1/nl. The first factor of n! comes from 
n! combinations which appear when the photon production operator act on the multiphoton 
state (k\, /c2-, k n \\(a^) n . The factor of 1/n! comes from the expansion of the action S = 
exp(i J d^xA^J^), and the factor of l/\fn\ comes from the normalization of the n-photon 
state. So the net result is proportional to 1/n!, which is exactly what is needed to get the 
Poisson distribution p n = exp(— n)n n /nl. In the case of the production of n-pairs, we have 
the same 1/n! from the expansion of the action, but now we get 1/n! coming from the 
normalization and not 1/y/nl as before. However, the action of the product of the creation 
operators of Q and L, which can be symbolically written as (ag6^) n , gives now an overall 
factor of nl from the action of, say, («g) n , as above, and also the sum of n! equal but 
not interfering terms, each of them being proportional to a different delta-function of the 
momenta, b(j>Q J + PL k )- Thus in the matrix element squared we will get the same overall 
factor 1/n! which is necessary for the Poisson distribution. 

Appendix C: Calculation of Baryon Asymmetry 

Here we calculate the lowest order nonzero contribution to the baryon asymmetry; we 
derive eq. (3.13) from eqs. (3.9) and (3.10). As our starting point, we have 




(C.l) 



where 




(C.2) 



and 




(C.3) 



Using eq. (3.8), we find that 



6(2u) 



9i r (-T/2 + im R ) 



(T/2 + im R ) 



(CA) 



Aiu L (~r/2 + im R + 2iu) 



(-r/2 - im R + 2ico) 
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and 



Thus 



02 f 2u,)* = r (imR + r/2) (-tm R + T/2) 

Aiu ^2im R + 2iu + T 2iu - 2im R + T 



2iu + T 



■] 



(C.5) 



16w 2 



{-ml - r 2 /4) 



+ 



(m 2 R - r 2 / 4 + Tim R ) 



+ 



(2im R + 2iu + T)(2iu + im R - T/2) (2iu - 2im R + T)(2iu + im R - T/2) 
T(im R - T/2) -m\ + iTm R + T 2 /4 



(2iu + T)(2iu + im R - T/2) (2im R + 2iu + T)(2iu - im R - T/2) 



m 



R 



+ r 2 /4 



T(im R + T/2) 



(2iuj - 2im R + T)(2iu - im R - T/2) {2iuj + T)(2iu - im R - T/2) 

(C.6) 

Now we must integrate each of the terms in eq. (C.6) as indicated in eq. (C.l). The lower 
limit of the integral is tuq + <C m R and we use T <C tuq + m^. We find that the 
first term cancels with its hermitian conjugate, the third and sixth terms are 0, the second 
and fourth terms cancel each other, and the fifth term plus its hermitian conjugate is 
responsible for the final result given in eq. (3.13), 



nB = n h -ni = 



16tt 



m R f 2 ef 



2/j3 



{CI) 



Appendix D: The Effects of Mixing in the Q, L Basis 

We will consider the decay of 6 to a QL pair (superscript 1 for this decay channel), 
and the decay of 9 to a QL pair (superscript 2 for this decay channel). For the first decay 
channel, from eq. (3.16) we see that a Q produced at the time t = is given by 

V>(0) = Q = sVl + cip 2 , (D.la) 

where 

' and s = - 6 = . (D.lb) 



Similarly, 

x(0) = Z = q?i-s^ 2 - (^-2) 

We will let the fields ip and x evolve in time, mixing their Q and L components as they 
travel. The time evolution of ip(t) can be modeled as follows: 

</,(£) = {se-^^i + cip 2 ) exp(-iuj 2 t) , (D.S) 
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where Auj = uj\ — U2- We now wish to ask the question: what is the Q content at some 
time t of the field ip which was initially pure Q? Using eq. (3.16), we can write eq. (D.3) 

as 

^{t) = [(c 2 + s 2 e- lAujt )Q - sc(l - e- lAut )L] exp(-tu 2 t) . (DA) 

The quark content is given by the magnitude squared of the coefficient of the first term, 
so that 

n ^\t) = [ c 4 + S A + 2c 2 s 2 cosAut]^ ^ / dk Q e% \A qL \ 2 . (D.5) 

sq,sz 

Similarly, from the same decay process 6 — > Q + L, the L that is produced can convert to 
a Q so that we have 



(t) = yJ2 2s2f;2 ( 1 " cosA ^) / dk Q dk L \ A QL\ 2 • ( D - 6 ) 



n ( - 1} 

SQ,Sl 



From 9 — > LQ, one can obtain Q at a later time from oscillations of either the L or the Q 
and find contributions: 

n {2 \t) = [c 4 + s 4 + 2c 2 s 2 cosAut]y ^ / ^ dhT L dkg \A l q\ 2 (D.7) 

and 



(2) 



(t) = Yl 2s2c2 ( 1 - cosAwt) y rf/cg |A X q| 2 . (D.8) 



s l,sq 

Thus the baryon asymmetry at any time t is 

n B (t) = n$\t) + n%\t) - ng>(f) - nf(t) 



[(c 2 - s 2 ) 2 + 4 S 2 c 2 cosA^] (\ A QL\ 2 ~ \A LQ \ 2 ) 

sl,s q (D.9) 



+ 4s 2 c 2 cosAcjt 



One can see that the baryon asymmetry oscillates in time as a cosine about the average 
value. When one takes a time average, the cosine term averages to zero, and one reproduces 
the result in eq. (3.28), 

Our derivation above assumes in eqs. (D.5-D.8) that all QL pairs and all LQ pairs were 
produced at the same time. If one considers that all pairs are not produced at the same 
time then an average over all pairs would also cancel the cosAcut term in eqs. (D.5-D.8). 
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